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Abstract 

A formulation of the Green's function method is presented in the n-pole 

approximation. Without referring to a specific model we give a general scheme 

of calculations that easily permits the computation of the "single-particle" 

Green's function in terms of the energy matrix. A theorem is proved which 

states that the moments of the spectral density function are conserved up to 

the order 2(n — Z + l), where I is the order of the composite field. A comparison 

with the spectral density approach is also discussed. 
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The theoretical analysis of highly correlated electron systems contains many unsolved 
problems and a complete and systematic formulation is still lacking. Many analytical meth- 
ods have been developed. Among many we recall the slave boson method [1], the non-crossing 
approximation [2], the doo method [3], the projection operator method [4], the method of 
equation of motion [5], the spectral density approach [6], the coherent potential approxima- 
tion [7], the composite operator method [8]. It is not easy to judge the reliability of the 
various approximations. Certainly one important aspect is the capability to conserve the 
symmetries inherent to the model and some relations, which can be derived on general basis. 

A scheme of approximation, common to various methods [4-6, 8], is the so called n-pole 
approximation, based on a linearization of the equations of motion. In this article we show 
that in this scheme the exact relation between the spectral moments and the spectral density 
is conserved: (i) ad, any order when the moments are calculated with respect to the linearized 
equations: (it) up to the order 2{n — I + 1), where / is the order of the composite field, when 
calculated with respect to the full equations of motion. This aspect is attractive and shows 
that the approximation can reproduce some essential features of the single-particle function. 

Let us consider a certain Hamiltonian 

H = H(ct> 1 ,...A P ) (1) 

the set {<f)i} denotes field operators. In order to close the infinite hierarchy of equations 
of motion some truncation is necessary. One procedure is to choose a basis of operators 
{■01, . . . ,ip n } and linearize the equations of motion as 

d - 

i— if>(x) = [i>(x), H] = e(-iV)^(x) + 6j(x) « e(-iV)ip(x) (2) 

where the eigenvalue or energy matrix e [9] is self-consistently calculated by means of the 
equation 

e(-iV x ){[^,t),^(y,t)]j = ([[^,t),H]^(y,t)]j (3) 

The symbol [. . .] ± denotes equal-time anticommutator or commutator, in dependence of 
the statistics of the set {ipi}. The rank of the energy matrix is equal to n, the number of 



components of the vector ip(x). When the basis {ipi} coincides with the original set {fa}, 
the linearization procedure given in Eq. (2) corresponds to the Hartree-Fock approximation. 
The fields ipi ma Y De called composite fields [8], since they are generally constructed as a 
combination of the original fields fa. 

Let us consider the thermal retarded Green's function [10] 

S(x,y) = (R [tP(x)^(y)}) = J d d kdwe^^-^-^ SQz,u) (4) 

where R is the usual retarded operator and the symbol (...) denotes the thermal average. 
d is the dimensionality of the system. By means of the Heisenberg equation (2), in the 
linearized form, the Fourier transform S^k, a;) is given by 

1 1 n_1 

where 

7(k) = F.r.([^(x,0,^ t (y,0] ± > (6) 

the symbol F.T. denoting the Fourier transform, and 

n 

D(k,tu) = J2^a t (k) (7) 

i=0 

is the characteristic polynomial of the matrix e(k). The characteristic coefficients Oj(k) are 
defined by the following relation 

a„_ fe (k) = (-l) fc Tr fc [e(k)] for < k < n (8) 

Where Tr k [e] is the trace of the fc'th order, defined as the sum of the determinants of all 
(j^j matrices of order (k x k) that can be formed by intersecting any k rows of e with the 
same k columns. We note that Tr n [e] = Det[e] and the convention Tr [e] = 1 is used. The 
A^(k) are the n x n matrices, determined by the equation 

iyAW(k) = D(k,u) [co - e(k)]- 1 J(k) (9) 

i=0 



The solution of this equation gives: 

n 

A(™)(k) = Yl a s (k)e s - m - 1 (k)/(k) for < m < n - 1 (10) 

s=m+l 

The Green's function (5) has a n-pole structure 

The energy spectra Ei(k) are the characteristic values of the matrix e(k) and correspond to 
the zeros of the characteristic polynomial 

n 

E Om(k)£T(k) = (12) 

m=0 

The spectral functions are given by 

i n— 1 

^) = rTu E W W (k) (13) 

°i{ K ) m =0 



where we put 

n 

Mk)= II [^(k)-^.(k)] (14) 
We also note that in the n-pole expansion, the Green's function can be expressed as 

to - £ Q/3 (k,u;) 

where the self energy £ a/ g(k, uj) is given by 

i n— 1 

Ectf(k,a;) = — Eu/[e(k)A»(k)U (16) 

E^g(k) 1=0 

i=0 

Summarizing, in this scheme of approximation the calculation of the Green's function 
requires the following steps: 

(%) Given a certain Hamiltonian, a basic set {■?/>«} of field operators is chosen; 
(ii) the energy matrix e(k) is calculated by means of Eq. (3); 

(Hi) the energy spectra Ei(k) are calculated as the roots of the characteristic polynomial of 
e(k); 



(iv) the AW(k) matrices are calculated in terms of the energy matrix e(k) by means of Eq. 
(10). 

We want now to study if the n-pole approximation conserves the exact relation between 
the spectral moments and the spectral density, that can be derived on general basis. We 
start by the following definition for the spectral moments 



M«(k) = F.T. 



id 
dU 



id 
dt„ 



ijj(x),^(y) ) 



[forO<p<A;] (17) 



We note that M^(k) is a n x n matrix and is a generalization of the quantity usually 
introduced in the literature, where the spectral moments are referred to the one-particle 
propagator. Let us introduce the spectral density function A(k, u) 



and consider the quantity 



A(k, u) = --ImS(k, u + irj) 



r+oo 

B {k \k) = / dujuj k A(k,uj) 

J — CO 



(18) 



(19) 



By taking time derivatives of Eq. (4) we find the following relation between the spectral 
moments and the spectral density 



M (fe) (k) = B (k) (k) 



(20) 



This relation is an exact one and must hold for any integer non negative value of k. By 
means of the field equation Eq. (2) the spectral moments can be written as 



M^ k \k) = M^ k \k) + SM^{k) 



(21) 



where 



M (fc) (k) = e fc (k)/(k) 



(22) 



<JAf(°)(k) = <5M«(k) = 



(23) 



5 



k-1 



SM (k) (k) = F - T - 



m=l 



,fc— m— 1 1 



id/dt x ) m {\6j(x),^(y)]) 



[for k > 2] (24) 



What see that the M^(k) correspond to the spectral moments calculated in the n-pole 
approximation, where the higher order field 5j(x) in the Heisenberg equation (2) is neglected. 
We will now prove that the spectral moments (k) satisfy the relation 

r+oc 



r+oc 

M^(k) = dutu k A(k,u) for any k 

J — oo 



(25) 



The proof is divided in two steps: (i) M^(k) and B^ik) satisfy for k > n the same 
recurrence relation; (ii) the equality M^(k) = B^(k) is satisfied for k < n. 

By recalling that a s (k) are the coefficients of the characteristic polynomial of the matrix 
e(k), the Hamilton-Cayley theorem says 



]Ta s (k)e s (k) = 
Then, the moments M( fc )(k) satisfy the recurrence relation 

n-l 

M n+k (k) = -J2 a s (k)M s+k (k) for k > 



(26) 



(27) 



s=0 



On the other hand, by noting that from Eq. (11) the spectral density is given by A(k, uS) 

n 

£)<7«(k)<S [u - Ei(k)], we have 



i=i 



i=i 



5( fc )(k)=^:^(k) ( T«(k) 

Since the E^k) satisfy the relation 

n-l 

E? +k (k) =-J2 a>m(k)E™ +k (k) for k > 

m=0 

the quantities B^(k) satisfy the recurrence relation 

n-l 

B (n+k) (k) = - a m (k)B {m+k \k) for k > 



(28) 



(29) 



(30) 



m=0 



We see that the spectral moments M^'(k) and the quantities B^(k) satisfy the same 
recurrence relation for k > n. To show that they coincide, it is enough to consider the case 
< k < n. 
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Let us take 1 < k < n — 1. By means of the expression (13) for the spectral functions 
we can write 



n-l 



B (k \k) = ]T \ (m) (k)P (k+m \k) (31) 



m=0 

where we have defined 

pt, ^-tM (32) 

By means of the relation (29) the quantities p( fc )(k) satisfy the recurrence relation 

n 

Y, a m (k)P m+k (k) = for k > (33) 

By simple algebraic relations it is easy to show that 

f0 for 0<fc<n-2 
P (fe) (k) = I (34) 
^ 1 for k — n — 1 

By recalling that from Eq. (10) A*™) = eA^ m+1 ^ + a m+1 I we can then write 

n— 2 

B (fc) (k)=A (n - 1) (k)P (Hn - 1) (k)+e(k) £ A (m+1) (k)P (fc+m) (k) (35) 

m=n— 1— fe 

n-2 

+/(k) 2 « m+ i(k)P (fe+m) (k) (36) 

m=n— 1— fe 

On the other hand A^ n_1 ^(k) = a n I(k), therefore the previous expression takes the form 

n-l 



P (fe) (k) =e(k) J B( fe - 1 )(k) + /(k) J2 a m +i(k)P (fc+m) (k) =e(k)p( fe - 1 )(k) (37) 

m=n— 1— k 

where the relation (33) for k > 1 has been used. But for k = it is easy to see that 
P(°)(k) = J(k), therefore 

P (fc) (k) = e fe (k)/(k) = M (fc) (k) for < k < n - 1 (38) 

Because of the recurrence relations (27) and (30), the equality (38) holds also for k = n and 
the step (ii) is thus proved. Then, the equality AP fc )(k) = B^(k) is true for any k and the 
relation (25) is demonstrated. 



We now turn to the more general question if the n-pole approximation conserves the full 
spectral moments M^(k), calculated by the exact Heisenberg equations. In the context of 
the Hubbard model it is known that in the 2-pole approximation the first 4 spectral moments 
relative to the electronic field are conserved [6,11]. For higher order moments the situation 
is not known. 

It is always possible to choose the basic set as 

d 1+1 
ij-Mx) = M*)> H} = Y, Ji P (-iV)ip p (x) [for 1 < Z < n - 1] (39) 
at P =i 

For this choice we have that 

e lv — llP 

forl<Z,p<n-l (40) 

Sji = 

and the matrix e has the following structure 

f I < n-2 

e ip = for (41) 
< p> 1 + 2 

From the definition (24) and by use of the exact Heisenberg equation (38) we can derive the 
following recurrence relation 

5M%\k) = e lr (k)SM%- 2 \k)el p (k) for 1 < l,p < n - 1 (42) 

Then, by means of Eqs. (24) and (40) we find that the matrix 5M^(k) has the following 
property 

5Mll k) (k) = 5M^ k+1) (k) = for 1 < I, p < n - k (43) 

By recalling Eq. (25) we can then enunciate the following theorem. 

Given a set of fields {ipi, I — 1, . . . , n}, if the subset {tpu I — 1, ■ ■ ■ ,n — 1} satisfies linear 

i+i 

Heisenberg equations % (d/dt) tpi(x) = ^2"fip(— iV)tp p (x), then the first 2(n — l + l) spectral 

P =i 

moments for the fields ipi [1 < I < n — 1] are conserved [12]. 

The theorem proved in this article shows that the n-pole approximation is equivalent 
to the SDA [6], when the polar ansatz is considered. The following remark is worth to be 



mentioned. In the SDA the unknown quantities, the energy spectra Ei(k) and the spectral 
functions crW(k), are evaluated in terms of the spectral moments through the equation 



To calculate one matrix element S a/ 3(k, cu) of the Green's function one must solve a set of 
2n non-linear equations. When n is not small this procedure is not very convenient. Indeed, 
most of the calculations in the SDA have been restricted to the case of two poles. In the 
n-pole approximation we only need to find the roots of the characteristic polynomial of the 
energy matrix. Once the energy spectra Ei(k) are calculated, the n x n Green's function 
can be immediately computed, as we have previously shown. 

In conclusion, in this paper we have presented a general formulation based on the n- 
pole expansion for the calculation of the "single-particle" Green's function. The denomi- 
nation "single-particle" is used in a generalized sense, and refers to the propagator of the 
n-component composite field {ipi}. Without recurring to a specific model, we have presented 
a scheme of calculations which allows us to determine the Green's function once the energy 
matrix e(k) is known. This is the fundamental quantity and contains the dynamics, con- 
trolled by the Hamiltonian, the algebra, determined by the statistics of the original fields 
{(pi}, the boundary conditions, since it refers to expectation values. The conservation of the 
spectral moments has been analyzed. We have shown that there is an internal consistency in 
the n-pole approximation, in the sense that the exact relation between the spectral density 
function and the spectral moments, when calculated on the same ground, is preserved at 
any order. We have demonstrated a theorem that states that the full spectral moments, cal- 
culated on the basis of the exact Heisenberg equations, are conserved up to a certain order 
that depends on the number of poles. For example, in tha case of a single-band electron 
model, the spectral moments of the electronic field are conserved up to the 2nth order. The 
method is equivalent to the spectral density approach. One advantage of this formulation 
with respect to the SDA is that we only need to solve the characteristic equation of e(k). 



n 




(44) 



i=i 
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